Oblique Half-Solitons and their Generation in Exciton-Polariton Condensates 
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We describe oblique half-solitons, a new type of topological defects in a two dimensional spinor 
Bose Einstein condensate. A realistic protocol based on the optical spin Hall effect is proposed 
toward their generation within an exciton-polariton system. 



Introduction. Exciton-polaritons (polaritons) are hy- 
brid quasi-particles that share the properties of both light 
and matter. On one hand, their excitonic part allows 
them to interact both with themselves and with the sur- 
rounding phonon bath. In microcavities, this peculiarity 
leads to efficient relaxation processes towards the ground 
state and to the formation of a Bose-Einstein conden- 
sate (BEC)[H with extended spatial coherence [2 ■ On 
the other hand, the photonic part of polaritons attributes 
them a very small effective mass which allows condensa- 
tion even at room temperature 

While the polariton BEC has indeed been obtained ex- 
perimentally, the superfluidity, namely the appearance of 
a dissipationlcss flow, was still an issue. Strong progress 
has been made in 2009 when the behavior of a flow of 
polaritons past an obstacle was analyzed 0, 01- One of 
the fascinating properties of a BEC is the possible cre- 
ation of topological defects by a perturbation. In spinor 
polariton condensates, elementary excitations are pecu- 
liar: indeed, two exciton spin projections (cr ± = ±1) 
can efficiently couple with light. Thus the wave function 
(order parameter) of a polariton BEC is no more scalar 
but vectorial. In such a framework elementary topologi- 
cal defects carry half-integer angular momentum and arc 
called half- vortices (HV). These objects were predicted 
theoretically^ in a thermodynamic limit, analyzed in 
presence of various (effective) magnetic fields 0-Q and 
observed experimentallvfioj MS M result disorder in mi- 
crocavities. 

A soliton is another example of topological defect 
that can basically be seen as a ID analogue of a vor- 
tex. It is characterized by a maximum phase shift of 
7r (dark soliton), on the scale of the healing length 
£ of the condensate, accompanied by a ditch in the 
density profile. In 2D, a new type of extended de- 
fect has been recently described: oblique dark solitons 
(ODS)[l3|, which were shown to be persistent or rather 
only convectively unstable [l2| ■ They can be generated 
if a supersonic condensate flow hits a defect larger than 
C[IH EH- These predictions have been recently verified 
in a polaritonic system fl4|. Polaritons can also exhibit 
bright solitons [15j. but these are not in the focus of the 
present Letter. On the other hand, solitons in spinor ID 
condensates 3, 13] and oblique solitons in spinor 2D sys- 
tems have also been theoretically considered 18j. In one 
dimension, many possible configurations were described, 



depending on the strength and type of the particle inter- 
actions (repulsive or attractive). In particular, a solution 
where the kink lies in only one component was reported: 
the dark-antidark soliton. This kind of defect reminds 
the HV configuration and we will see why it deserves to 
be called a " half-soliton" [Tg| . 

In the present work we discuss the generation of 
oblique dark half-solitons (ODHS), in a spinor BEC. The 
paper is organized as follows. First, we analyze the ID 
half-soliton by analogy with the half- vortex. Second, we 
introduce the ODHS. Third, we propose two experimen- 
tal configurations leading to their formation, capitalizing 
on the specificities of polariton condensates. The first 
protocol is using a defect acting only on one circular com- 
ponent. The second one is using the TE-TM splitting 
present in planar microcavities. 

The ID half-soliton. A two components polariton 
condensate at OK can be described by a spinor Gross- 
Pitaevskii equation [2pj . Assuming first a parabolic dis- 
persion with an effective mass m* and an infinite lifetime 
of the particles, one has 
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Time-independent solutions are found upon expressing 
the BEC wave functions as ^>± (r, t) = %p± (r) exp (—ifj£), 
where u = (a\ + a^^o is the chemical potential re- 
lated to the density of the homogeneous condensate 
no = no+/2 = "0-/2. Indeed, for a realistic polariton 
BEC, the interactions between the particles of opposite 
spin are weak and attractive: «2 — — 0.2a\ which leads 
to a linearly polarized ground state. When looking for a 
single- kink solution such as a vortex or a soliton, it is con- 
venient to discuss the asymptotic behavior. We assume 
that far away from the defect the condensate density is 
unperturbed. 

The order parameter can be written in two distinct 
representations 0, Q. On the basis of circular polar- 
izations discussed previously, each component possesses 
its own phase 9±(r) and (ip + ,if)-) = ^/n^ll (e l6l +,e ie -). 
In the linear polarization basis, the polarization angle 
rj(r) and the global phase of the condensate <p(r) can 
be separated: (ip x ,ip y ) = ^/uq (e IV cos (77) , e lip sin (77)) . 
The transformation from one to another is obtained via 
ip± = {ipx T tyx) /V% and 0± = <-p T V- I n the latter 
form, half-vortices are characterized by two half-integer 
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winding numbers for both the phase and the polarization 
which lead to integer shifts of 7r around their core. Similar 
considerations can be applied to solitons. Indeed, in one 
dimension the scalar dark soliton solution [2l[ is simply 
given by ips (%) = y/no tanh (x) and its phase is a Heavi- 
side function of amplitude n. In the circular polarization 
basis and, in the simplest case where ai — 0, a half- 
soliton (HS) is nothing but a usual soliton appearing in 
one component (let's say ip-) while the other remains ho- 
mogeneous. Thus, the associated order parameter reads 
(ip+ s ,ip^ s ) = y/rio/2 (1, tanh (a;)). Rewriting the latter 
on the linear polarization basis leads to (ipx S ,ipy S ) 



/uq/2 (1 + tanh (x) ,i — i tanh (x)). Looking at asymp- 
totic forms, one easily obtain 
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with h and s half-integer numbers which can be seen 
as topological charges. Basic HS appear for (h, s) = 
{(±1/2, ±1/2), (±l/2,=Fl/2)} and their phase and po- 
larization angle are ranging from to ir/2 with a fully 
circularly polarized center just like the half-vortex core. 
This topological defect can also be seen as a domain wall 
with respect to x- and y-polarized particles. A plot of 
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±,x,y) is 



the HS density profiles (rij 
proposed in FigQJa). 

The oblique half-soliton. For a 2D condensate it was 
shown 11. 13, 0] that stationary extended oblique dark 



solitons can occur both in the scalar and spinor sys- 
tems despite their instability, because the latter is only 
convective. They appear in a perturbed fluid at super- 
sonic velocities as merging vortices arranged in "trains" 
("streets"). 

Let us now focus on the possibility of creating ODHS. 
First of all, it is clear that in the case where the two 
components of the spinor BEC do not interact («2 = 0), 
if they are initially equally populated, a perturbation in 
only one of the components will lead to the formation of 
half-integer excitations. Next, what happens if the inter- 
action is no longer negligible? To answer this question 
we turn back to Eqs.(l) and (2), rescale them, look for 
stationary solutions and switch to the hydrodynamic pic- 
ture of Rcfflil where the phase of each components is ex- 
pressed by means of their stationary and irrotational ve- 
locity fields via v± (r) = h/m*V8± (r), with r = (x,y). 
We look for oblique solutions that depend only on the 
tilted coordinate x = (x — ay)/vl + a? which leads to 
the following set of equations 

(< 2 /4 - n CT </2) + 2n\ (Am, + 2A 2 n_ ff ) 

= (q + 2ii)nl -qnl (3) 
where Ai 2 = Oi\,7,/{oii + a-i) and q = U 2 / (1 + a 2 ) (U the 
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FIG. 1. (Color online) (a) HS density profiles scaled to no 
and £. The solid blue and red curves represent the n+ and 
n- density profiles while the dashed-dotted purple and dashed 
cyan curves show n x and n y respectively, (b) ODHS density 
profiles normal to its axis. The solid blue and red curves show 
numerical profiles, the dashed black curve is the perturbative 
solution described in the text and the dashed dotted green 
curve show the scalar soliton solution. The parameters are 
a = 10, U = 3.5, Ai = 1 and A 2 = -0.2. 



velocity of the flow). This system has to be solved nu- 
merically, but let us start with some simple arguments. 
The density profile of an integer oblique dark soliton [ll| 

is given by n DS = 1 - (1 - q/n) sech[xVM ~ l] witn 
/x = AiUq. Now, for the case of the ODHS, the den- 
sity notch in the o~ component, that contains the de- 
fect, is seen as an external potential by the initially un- 
perturbed er + component, because of the interactions 
between the particles of different spins. We suppose 
that the <r + component fits the shape of this potential 
which is nothing but ^n-. Then, this perturbation cre- 
ates in turn a potential for the a~ component given by 
— A2TI+ = —A 2 n-. Therefore, the density profile is mod- 
ified as n_ (Am- — AjTi-) /Ai. Iterating this proce- 
dure leads to a geometric series and to a renormalization 



of the interaction constant Ai <— Ai — Ai,/ {A\ — A2). 
Consequently, the ODHS solution is approximated by 



noDHS = 1 - (1 - <?/m) secn XVV - q 



(4) 



with /i = Aino- In this description, the sound velocity is 
changed like c s — > c s = y/ji/m* and the healing length 
like £ — > £ = h/ y/2m*Ji. In the case where A2 < 0(> 0), 
which corresponds to an attractive (repulsive) interac- 
tion, c s is slightly increased (decreased) and inversely for 
£. The component without a soliton obviously presents a 
minimum (maximum). This argumentation is compared 
to direct numerical solutions of Eqs. (4) and (5) in the 
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FigQJb) showing a remarkable accuracy for small values 
of A 2 . 

Generation of ODHS in a polariton fluid. Bringing 
an external perturbation in a quantum fluid is an effi- 
cient way to nucleate topological defects in a controlled 
manner. A relevant experimental setup relies on the ac- 
celeration of a superfluid above the sound velocity past 
an obstacle larger than £, where the breakdown of su- 
perfluidity allows the elementary excitations to develop. 
The defect is a potential barrier resulting from an off- 
resonance laser beam. A planar microcavity is a system 
that offers a high degree of control on the velocity and 
the density of 2D spinor BECs as discussed in Refs(3 
and El 

Now, to describe more accurately the polariton con- 
densate, we take into account the real non parabolic dis- 
persion of the particles via a set of four coupled spin- 
dependent equations for the photonic $f± and excitonic 
fields respectively. A resonant pumping situation is 
considered with photon lifetime of 100 ps. 
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The new quantities that appear are a quasi-resonant cw 
pumping which is taken to be linearly polarized P + = P_ , 
the strong coupling between exciton and photons with 
a Rabi splitting f2# = 10 meV and the photonic TE- 
TM splitting [23J of strength j3 for which we assume a 
quadratic dependence over the wave vector Q. We as- 
sume that the potential barrier we consider D± can affect 
independently each circularly polarized photonic compo- 
nent. 

In the numerical experiments that follow we repro- 
duce the protocol described in Refs|l3j and [lj involv- 
ing an immobile defect and a moving fluid. The pump 
laser is tuned slightly above the bare polariton branch 
(+0.3 meV) to remain resonant despite the interaction- 
induced blue shift. The spot (15 zzm upstream from the 
defect) is bar shaped (x x y = (10 x 100)/xm 2 ). The 
fluid propagates from the left to the right with a wave 
vector fc p ;r = 1/im -1 , low enough to avoid parametric 
instabilities, and collides against a 6 /zm large circular 
defect. Defining the sound velocity is far from trivial so 
far as the local pumping and propagation effects induce 
density fluctuations and its general exponential decrease. 
The sound velocity c s = sj ji/m* of an unperturbed flow 
is expected to decrease like -y/exp (—m*d/hkT p f l ) where 
77i* ~ 2771*^ and d is the distance from the pump. Strictly 
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FIG. 2. (a) Pseudospin representation of the order parameter 
at M ~ 1.8. The colormap shows the S z component (de- 
gree of circular polarization), (b) Density slice displaying 71+ 
(blue curve) where the soliton appears, n_ (red curve) with 
a minimum at the soliton position and the sum of the two 
(dashed black curve). The dashed white lines indicate the 
slice position. 



speaking, c s has to be defined locally and separately for 
each components. We choose a pseudospin representa- 
tion, which is relevant to provide information about both 
components of the order parameter in a single picture. 
We recall that the pseudospin S (r, t) is a 3D vector ly- 
ing on a Poincare sphere. Its components are defined 
by 



5ft (f+V- * 



Sy = S$ (iPf*^ 

S = ( '^P h - 



(7) 



normalized to unity. S x . y define linearly polarized states 
while S z is the degree of circular polarization of the 
condensate [23|. Below, we propose two different experi- 
mental configurations for the generation of ODHS. 

The first alternative is to find a way to perturb only 
one component. This can be done for polaritons where 
the defect potential could be optically created by a cir- 
cularly polarized pulse [2^. This scheme is however far 
from being ideal since it would bring of lot of unwanted 
perturbations to the system and it should rather be seen 
as a model experiment. In this framework, we impose 
D + = 20 meV, L>_ = 0, and a 2 = -0.2ai. The S z (x,y) 
component of the corresponding stationary flow regime, 
is presented in FigJ^a) for a Mach number M ~ 1.8. The 
generated pair of oblique solitons is left circularly (cr_) 
polarized and with the help of a density slice normal to 
the flow FigfJIh) we are able to undoubtedly identify an 
ODHS. 

Our second proposal is, on the other hand, completely 
realistic. The impenetrable defect is restored in both 
components and we suggest to benefit from the spatial 
separation brought by the k-dependent TE-TM splitting. 
The strength of the induced effective magnetic field £Ilt 
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is given by (3 = h 2 (mj~ — m^ 1 ) /4 where mi,t are the ef- 
fective masses of TM and TE polarized particles respec- 
tively. The energy splitting between the two latter eigen 
modes is Alt = Pk 2 x — 0.02 meV for n%t = 5 x 10~ 5 too 
and mi = 0.95m t (mo is the free electron mass). More- 
over, the field makes a double angle with respect to the 
direction of propagation and is thus symmetric with re- 
spect to the x-axis. The latter leads to the so-called 
optical spin Hall effect: the precession of the pseudospin 
around Q r t [25L l2(i| . In our numerical setup, the linearly 
x-polarizcd injected particles are initially moving along 
the ^-direction so no polarization separation is expected 
until the fluid reaches the obstacle. Indeed, at this pre- 
cise moment, the flow is split into two parts moving in 
oblique directions. These two parts start to precess about 
the effective field in opposite directions. Particles going 
to the left gain a a + component. Particles going to the 
right gain a er~ component. This antisymmetry results in 
the formation of an oblique soliton in one circular compo- 
nent, and not in the other, which is an ODHS. The figure 
[3](a) shows the resulting polarization pattern achieved in 
the steady state: A pair of oblique dark solitons of oppo- 
site circularity is formed behind the obstacle. It should 
also be noted that the Cerenkov radiations exhibit a re- 
markable polarization separation upstream from the de- 
fect. The figure [3Jb) shows a density slice of each circular 
component that reveals the ODHS pair. The profiles are 
qualitatively similar to the one obtained for the ideal case 
presented on the figure [Tfb) . The [3jc) shows the results 
of a similar calculation, but for a splitting twice larger. 
One can clearly see the formation and disappearance of a 
pair of ODHS close to the defect, and, further away, the 
formation of a second one of opposite circularity. The 
figure [3jd) shows the density when the four ODHS are 
simultaneously present. This behavior is due to the TE- 
TM splitting which couples the circular components and 
therefore couples ODHS of opposite circularity. To finish, 
we make the following remark: In the transient regime, 
before the steady state is established or at lower fluid ve- 
locities, we have clearly observed the generation of HVs, 
and it will be the topic of a future work. 



Conclusions. We have reported a new type of half- 
integer topological defect that occurs in two-components 
2D BEC: the oblique dark half-soliton. We have shown 
that an exciton-polariton fluid is well-suited for the study 
of such defects. We have made a realistic proposal, based 
on the optical spin Hall effect, which should lead to their 
experimental observation. Besides, the setup we suggest 
should also allow a controlled nucleation of half- vortices, 
which was never reported so far. 
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FIG. 3. Same representations as Fig[2] ODHS pair(s) in the 
presence of the TE-TM splitting, (a)-(b) f3k^ x = 0.02 meV 
and (c)-(d) fik 2 v . x = 0.04 meV. 
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